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Abstract 

We investigate the single-impurity Anderson model by means of the recently introduced modified perturbation theory. 
This approximation scheme yields reasonable results away from the symmetric case. The agreement with exactly known 
results for the symmetric case is checked, and results for the non-symmetric case are presented. With decreasing 
conduction band occupation, the breakdown of the screening of the local moment is observed. In the crossover regime 
between Kondo limit and mixed- valence regime, an enhanced zero-temperature susceptibility is found. 
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1 Introduction 

The single-impurity Anderson model (SIAM) was 
originally introduced to describe magnetic impuri- 
ties in weakly correlated non-magnetic metals JjJ. It 
was successfully applied to explain unexpected low- 
temperature thermodynamic properties and photoe- 
mission data ^ [| . Although its predictions for pho- 
toemission in periodic structures is still a matter of 
controversial discussions |5|, ^j], it remains never- 
theless surprisingly well suited for describing alloys 
of magnetic and non-magnetic compounds. 

The SIAM is one of the most-examined and best- 
understood models of many-body physics. The renor- 
malization group theory and the Bethe-Ansatz 
(see e.g. ||) have contributed enormously to the un- 
derstanding of its physics. An excellent review of 
these and other theoretical approaches to the SIAM 
can be found in ||. 

In the recent past, the SIAM has gained much in- 
terest in the context of the dynamical mean-field the- 
ory (DMFT) 
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It has been shown that in the 
limit of infinite dimensions, correlated lattice models 
can be mapped onto single-impurity models like the 
SIAM. Since the Bethe ansatz method requires spe- 
cial assumptions concerning the conduction band, it 



is not applicable in the DMFT. Often, numerically 
exact methods such as exact diagonalization (ED) or 
quantum Monte Carlo (QMC) are used. Besides of 
being numerically expensive, there are also physical 
limitations to these methods. ED is principally lim- 
ited to finite sized systems, while using QMC it is 
difficult to obtain results for low temperatures. Both 
methods are not well suited for calculating spectral 
densities. In addition to these methods, numerical 
renormalization group calculations to both the SIAM 
and via DMFT to the Hubbard model have been pre- 
sented This method gives by construction very 
useful information for the low-energy behaviour, but 
the quality concerning high-energy features remains 
an open question. So there is still a need for reliable 
analytical methods. 

A standard approximation scheme is second-order 
perturbation theory around the Hartree-Fock solu- 
tion (SOPT-HF) @. In the context of the DMFT 
it is called iterative perturbation theory (IPT) |TT) . 
This method gives convincing results at special pa- 
rameter sets for the lattice problem, which map onto 
a particle-hole symmetric SIAM, but does not work 
very well for arbitrary system parameters. To extend 
the IPT to arbitrary band-fillings, an ansatz based on 
the SOPT-HF self-energy was introduced 111 pi Il6), 
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which reproduces the atomic limit and the Friedel 
sum rule [jl7|, [l8[ [l9| as well. On the basis of 
this ansatz, a further improvement was recently pre- 
sented |2l[] using the first four moments of the 
spectral density to determine the free parameters. 
The method (modified perturbation theory MPT) 
works away from the symmetric point and the mag- 
netic properties of the solution |22| agree well with 
QMC results ^3|. Since these results were obtained 
for the Hubbard model within the DMFT, it is still of 
interest to solve the original SIAM within the MPT. 
After outlining the main ideas of the MPT in the next 
section, we will discuss the results for the SIAM both 
in the symmetric and non-symmetric case. Compar- 
ison with exactly known results will provide a good 
basis to estimate the quality of the MPT. 

2 Theory 

The Hamiltonian of the SIAM reads as 



Cko 



k,a 



(1) 
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where cj}^ annihilates (creates) a conduction band 
electron with quantum numbers k and spin a, e& is 
the corresponding eigenvalue of the conduction band 
Hamiltonian and /i the chemical potential, c d j^ anni- 
hilates (creates) an electron at the impurity-level e d . 
"■do- = c da c d a is the impurity occupation number op- 
erator. The hybridization strength between conduc- 
tion band and impurity electrons is given by Vk d and 
enters all practical calculations via the hybridization 
function 



A(£) = ]T 



V 2 
v kd 



E + [i - e k 



V 2 dx 



Po{x) 
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In this paper we assume the hybridization strength to 
be constant (Vkd = V), although this is not required 
for the MPT. po(E) is the free density of states of 
the conduction electrons with bandwidth W. The 
impurity electrons interact with each other, where the 
interaction strength is U. Contrary to lattice models, 
e. g. the Hubbard model and the periodic Anderson 
model, the interaction takes place only at one single 
impurity level. 

All relevant information can be obtained from the 
impurity Green function G da (E) = ((cd CT ; e f° r 



which we find the formal solution using the equation 
of motion 



G drT (E) 
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E + (i - e d - A(E) - £ dCT (#) 



(3) 



where we introduced the impurity self-energy ~E da (E) . 
The actual problem is to find an (approximate) so- 
lution for the self-energy. We choose the modified 
perturbation theory [p0| , which was applied to the 
Hubbard model pi[ |22| within the framework of the 
dynamical mean- field theory fl(i| , 11 1. Starting point 
is the following ansatz for the self-energy |Q : 



£ d(T (£) = U(n d - a ) + 
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where S^ 00 - 1 (E) denotes the second-order contribu- 
tion to the SOPT-HF. The parameters a a and b a are 
chosen such that the first four spectral moments (up 
to n = 3) 
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are correctly reproduced. 

A da (E) = (-l/7r)3G dc r(£ + iQ + ) is the spectral den- 
sity. This requires pOfl: 



(n d _ g )(l - (nd-o-)) 
(n d _ CT )( HF )(l-(n d -a)( HF )) 
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where the superscript (HF) denotes Hartree-Fock ex- 
pectation values. B dlJ is a higher correlation func- 
tion introduced by the fourth moment, which can be 
expressed in terms of dynamic one-particle quanti- 
ties m 

{n da ){l - (n d a))(B dc r - e d ) = 

= E^(cLcd,(2n d - CT -l)) (8) 

k 

= ~S J dE f_(E)A(E) ^da(E) - lj G da (E) 

with the Fermi function f-(E) = (exp(pE) + l) -1 . 
The results from reference p6| are simply recovered 
by replacing B da by its HF value: 

(^} (HF) (i-w (Hr) )(5r ) -^)= 

= E^( c L^) (HF) (2(n d - CT ) HF -l)) (9) 

k 

= ~Qj dEf_(E)A(E) (2(n d - CT ) (HF ) - l) G d f ] (E) 
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The fictitious chemical potential pt CT that appears in 
equation (0) is introduced as an additional parameter 
in the HF Green function 



E + ri a -e d -A(E)-U(n d - a ) 



(10) 



It can be used to enforce for instance the Friedel sum 
rule as was done in jl6j , but one would restrict oneself 
to T = 0. Therefore, we choose the condition (n dcr } = 
(«do-)' HF ' ) to fix Jl a pC| . Thus, (0) is always simplified 
to a a = 1. 

The modified perturbation theory constructed in 
that way reproduces several non-trivial limiting 
cases || such as U/W^ -> (W^ = V 2 /W 
is the effective bandwidth of the impurity-level due 
to the hybridization) and W^* 1 /U — » (in the sense 
of H). 

In the symmetric case (i. e. symmetric density of 
states, half filled conduction band and impurity level, 
e d = —U/2) b a vanishes, and the ansatz (||) re- 
duces to the SOPT-HF, since the latter reproduces 
the first four moments accidently. It is known, that 
the SOPT-HF in the symmetric case gives reasonable 
results even for large U ]25| . 

3 Results and discussion 

3.1 The symmetric case 

First we like to summarize the results for the sym- 
metric case. In this case, the chemical potential will 
always be fixed at the center of gravity of the con- 
duction band. In our calculations, the conduction 
band density of states is semi-elliptic and of unit 
width centered at E = 0, thus defining the energy 
scale used throughout this paper. In figure [l] (upper 
picture) , the impurity quasiparticle density of states 
(QDOS) is plotted for various values of the hybridiza- 
tion strength V with constant [7 = 1. For all values 
of V, it consists of three structures. Two are located 
approximately at — y and y, respectively. These cor- 
respond to the atomic quasiparticle levels (e d , e d + U) 
and are commonly called charge excitations. The cen- 
tral peak positioned at /i is usually denoted Kondo 
resonance (KR) |l. This peak is a genuine many- 
body effect. The states responsible for this feature in 
the QDOS are ascribed to antiferromagnetic correla- 
tions between conduction electrons and the localized 
electrons on the impurity site which form a local mo- 
ment. The antiferromagnetic correlations provide a 
screening of the local moment (Kondo screening) |^| . 
Both the width 5k and height of the KR scale with 
V. The inset in figure [l] shows the variation of 5k 
with V. 




Figure 1: Impurity quasiparticle densities of state 
(impurity QDOS) of the symmetric SIAM (e d = 
n d = 1.0, n c = 1.0, T = 0). Top: U = 1.0 and differ- 
ent hybridization strengths V; Bottom: V = 0.2 and 
various U. Only the vicinity of E = \i = is plotted. 
Inset: width 5k of the Kondo resonance. 



The Kondo resonance is plotted in the lower part of 
figure [l] for different values of U and constant V. The 
height is constant, since for the symmetric case, the 
MPT fulfills the Friedel sum rule as will be discussed 
below. The width 5k varies as ~ 1 _ co 1 nstu -2 ■ 

The impurity susceptibility ^( im W (T) is defined as 
the change of the susceptibility at the impurity site 
(denoted by i) due to the presence of the impurity: 



X 



(imp) ^ji^ _ 



<9((n dT + n ci] ) - (n dl + n ciL )) 



dB 



B-»0 



X 



(11) 

(o) (T) 



with xf^ being the susceptibility per lattice site of 
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Figure 2: T X {lmp) (T) and x (lmp) {T) (inset) as func- 
tions of the temperature for the symmetric SIAM 
for different values of U (V = 0.2). The x-axis of 
the T X ^ mp \T) -plot is normalized to the respective 
Kondo temperature (see text). 



the pure system, i.e. without the impurity, and 
n cla = c\ a c la . Figure | shows T X {lmp) {T) as a func- 
tion of T. The temperature scale is normalized to the 
Kondo temperature defined by Tk = ^p^5T7t=o) © ' 
This quantity is proportional to the inverse width 
of the KR: T K ~ ^ @. For large values of U, a 
scaling behaviour is observable, the susceptibility as 
function of is approximately independent of U. 

For high temperatures, x ^ lmp ^ (T) becomes Curie-like: 
X (imp) (T) - i. Thus the impurity behaves like a lo- 
cal moment of size 

Tx (im P ) ( T ) _ For low temperatures, 
strong deviations from the Curie law are obvious. For 
small U : x ( im P)(T) becomes flat, showing Pauli-like 
behaviour. For U > 2.0 (Kondo limit), the suscep- 
tibility develops a maximum, the position of which 
scales with Tk and marks the onset of the screening 
of the local moment. For T — > 0, the susceptibility is 
Pauli-like again, but with reduced value. 

In summary, for the symmetric SIAM the MPT cal- 
culations, which are equivalent to SOPT-HF in this 
case, agree well with expected or exactly known re- 
sults ||. This is of course due to the fact, that for 
the symmetric SIAM, the perturbation series for the 
self-energy converges rapidly even for large U pa]. 



3.2 The asymmetric case 

In this section, we present results for the non- 
symmetric SIAM. Away from the symmetric point, 
less results are known exactly. Furthermore, stan- 
dard perturbation theory cannot give as convincing 
results as in the symmetric case. 
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Figure 3: Impurity QDOS of the non-symmetric 
SIAM: U = 2, e d = -1., V = 0.2, T = 0, but various 
n c . The chemical potential is positioned at the ar- 
rows. The inset shows the chemical potential \x and 
the susceptibility y^ lm P)(T — 0) as function of the 
conduction band occupation n c . 



First we study the variation of the conduction band 
occupation number n c . Starting from the symmetric 
case with U = 2 and V = 0.2, we vary n c between 1 
and 0. Correspondingly, the chemical potential shifts 
towards lower energies. In figure ||, the QDOS for 
different n c is plotted. The arrows indicate the po- 
sition of the chemical potential /i. The inset shows 



X 



(imp) frp 



(T = 0) and \i as function of n c . With de- 
creasing n c , the Kondo resonance also shifts to lower 
energies. Note that for n c < 1 the chemical potential 
lies below the center of the KR, which additionally 
becomes rather asymmetric. For very low carrier con- 
centration n c , the KR disappears. For most carrier 
concentrations is the zero-temperature susceptibility 
almost constant, only a slight increase with reduc- 
ing n c is visible. At n c w 0.07, however, ^( irn W(0) 
diverges. This divergence is directly related to the 
disappearance of the Kondo resonance. Exactly at 
n c = 0.07, the chemical potential /i crosses the lower 
edge of the conduction band. There are no electrons 
which could form (via antiferromagnetic correlations) 
a Kondo resonance, thus no screening occurs. The 
still occupied impurity states of the lower charge ex- 
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Figure 4: Impurity QDOS for different values of ed 
and constant U = 4, V = 0.2, n c = 1.0 (therefore 
tx = 0). 



citation peak form a local moment. Since now ^("™p) 
is the susceptibility of a stable local moment, the di- 
vergence is to be expected. Only for n c = the sus- 
ceptibility vanishes completely. This is, because the 
chemical potential shifts to fi — » — oo, and the impu- 
rity level is therefore also empty. 

Next, we like to investigate another asymmetric sit- 
uation: we fix U to a large constant value of U = 4 
and vary the energetic position of the impurity ed 
with respect to the conduction band. The hybridiza- 
tion strength is set to V = 0.2 and the occupation 
number of the conduction band n c = 1. The latter 
implies fi = 0. The situation is visualized by means 
of the impurity QDOS in figure ||. The lower charge 
excitation peak moves with increasing ed from below 
(ed = —0.75) into the conduction band. The upper 
charge excitation at approximately ed + U also moves 
towards higher energies. Contrary to the charge ex- 
citation peaks, the Kondo resonance remains at its 
position E = fj, = 0. For ed > the KR merges into 
the lower charge excitation peak. The upper charge 
excitation looses spectral weight as the impurity level 
becomes depopulated. For ed > 0.5, a dip indicates 
the onset of a gap within the lower charge excitation. 
This gap will simply separate the charge excitation 
from hybridization-induced states within the conduc- 
tion band and is not a correlation effect. 
In particular, figure 4 shows that within the MPT 
the charge excitation peaks are positioned close to 
the atomic quasiparticle levels as expected. Con- 
trary, within the SOPT-HF the positions do not sig- 
nificantly depend on ed- 
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Figure 5: ^(""^(T) for the parameter sets of figure 
I 



The transition from the Kondo limit to the in- 
termediate valence regime (IV) and finally to the 
empty-impurity state manifests itself in the suscep- 
tibility (see figure^). For high temperatures the sus- 
ceptibility follows Curie's law, X (mp) [T) ~ ? in all 
three cases. Only for low temperatures, y^ lm P)(T) be- 
haves differently for the respective parameter regimes. 
In the Kondo limit (ed < —0.5, solid and dashed 
line)), x ("™p)(T) shows a maximum and decreases 
for decreasing temperature. For very low T < 
O.ITk an increase is visible, for even lower tem- 
peratures, the susceptibility becomes constant again 
and ^x < -™ p H'^)|t=o = (not recognizable in figure 
5). In the IV-regime, and most pronounced around 
ed ~ —0.25, the y-behaviour persists down to much 
lower temperatures, the zero-temperature susceptibil- 
ity is enlarged, and therefore the Kondo temperature 
Tk becomes smaller. This can be ascribed to a weaker 
screening of the local moment. A further increase 
of ed depopulates the impurity. The "effective local 
moment" shrinks and the susceptibility gets reduced. 
The T-axis in figure |^ is normalized to Tk using the 
definition mentioned in section 3.1. No scaling be- 
haviour is visible. 

The above-mentioned enhancement of ^(™p) (T) in 
the Kondo regime for low temperatures, which is most 
clearly observable for ed = —0.5 needs further discus- 
sion. This feature is obviously most pronounced at 
the crossover from the Kondo limit to the IV-regime. 
An interesting observation can be made by examining 
the impurity density of states pd(E) at E = fi = 0. 
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Figure 6: Lower part of the impurity QDOS for 
e d = -0.5, U = 4, V = 0.2 and n c = 1.0 at vari- 
ous temperatures. The inset shows the value of the 
impurity QDOS p d (E) at E = p, = 0. 



For ed = —0.5 (see figure 6), it shows a maximum 
at the same temperature where the increase of the 
susceptibility begins. For e d away from —0.5, this 
maximum becomes less pronounced and eventually 
disappears. With its disappearance the characteris- 
tic increase of the susceptibility also vanishes. This 
as well as the concurrence of the position of the max- 
imum and the onset of the increase of the suscep- 
tibility suggest a closer relation of the behaviour of 
Pd{p) and x^ mp \T). However, it is astounding that 
a decrease of /9 d (p) provokes an increase in ^( im P) (T) . 
One possibility to explain this behaviour is the follow- 
ing: The states, that build up the Kondo resonance, 
are responsible for the screening of the local moment. 
In the examined non-symmetric case, most spectral 
weight of the KR lies above p (see figure (||)). Only 
those electrons filling the states close to the chemical 
potential contribute to the screening. If Pd(p) goes 
down, the number of electrons participating in the 
screening decreases. The screening should become 
less effective. Eventually, the local moment, which 
is still present, shows a stronger reaction to an ap- 
plied magnetic field. Consequently, the susceptibility 
increases. In this way, the observation may be under- 
stood. 

To finish our survey of the asymmetric SI AM, we 
like to check the Friedel sum rule B [lTj numerically. 
The Friedel sum rule relates the occupation number 
of the impurity level, (ndo-), to the real part of the 
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Figure 7: Test of the Friedel sum rule for the asym- 
metric SIAM. The impurity occupation is plotted as 
dotted line, the result for the impurity occupation 
using the Friedel sum rule as solid line (dashed line: 
SOPT-HF result). Left picture: U = 1.0, V = 0.2, 
e d = -0.5; right picture: U = 1.0, V = 0.2, n c = 1.0. 



self-energy at the Fermi energy: 

(n.) 1 - i arctan f fi + ^l + ^O) \ 
(n dCT )- 2 ^ arctan^ J 

(12) 

This formula can be derived using the Luttinger the- 
orem Q as shown in || pl|. 

For the symmetric SIAM, (Q) becomes trivially 
fulfilled in the MPT (equivalently in the second order 
perturbation theory) since 9£A(0) = 0, 5RS CT (0) = f 
and e d = — -j. Therefore the impurity occupation is 
always (ri d(T ) = 0.5. Our results for the asymmet- 
ric case are plotted in figure |?]. The dotted line is 
the impurity occupation (n dcr ) within our approxi- 
mation. The solid line (dashed line) is obtained us- 
ing @ with the MPT (SOPT-HF) self-energy. In 
the left picture, the conduction band occupation n c 
is varied from 0.3 to 1. The other parameters are 
chosen such that the right boundary represents the 
symmetric case for {7 = 1. As discussed above, in 
the symmetric point both the MPT and the SOPT- 
HF fulfill the Friedel sum rule. For n c < 1, however, 
the SOPT-HF deviates almost immediately, whereas 
in the MPT, the Friedel sum rule is met in excellent 
agreement until n c « 0.55. For lower n c , the devia- 
tion increases. The origin of this breakdown within 
our theory is not yet clarified. The right picture shows 
the situation of fixed values of U = 1 and n c = 1, e d 
being varied. For the left boundary, the parameters 
are identical to the right boundary of the left picture 
(symmetric point). Varying e d away from the sym- 
metric point both the MPT and the SOPT-HF do not 
fulfill the Friedel sum rule. However, in this case, the 
superiority of the MPT relative to the SOPT-HF still 
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manifests itself in the correct positions of the charge 
excitations (figure 4). The SOPT-HF fails to repro- 
duce these. 



4 Conclusions 

We have analysed the SIAM within the modified per- 
turbation theory (MPT). For the symmetric SIAM, 
the MPT and therefore our results are equivalent to 
the SOPT-HF. In this case, the SOPT-HF gives rea- 
sonable results even for large values of U . For the 
asymmetric case, the SOPT-HF does not work very 
well. This may be due to the fact that the SOPT-HF 
fails to reproduce the first four moments of the spec- 
tral density. However, the MPT, which by construc- 
tion always fulfills these moments, gives reasonable 
results also away from the symmetric point. Within 
the MPT, the positions and weights of the charge ex- 
citations are predicted correctly since this is closely 
related to the first four moments p^ j. Furthermore, 
we achieve a better agreement with the Friedel sum 
rule compared to the SOPT-HF. 

Concluding, one can state that the MPT yields 
reasonable results in a wide region of the parame- 
ter space. Its conceptual simplicity and its numerical 
feasibility qualify the MPT as a valuable complement 
to the QMC method in the context of DMFT. 
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